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We show that a composite spin-y 2 Kondo system can be formed by coupling a superconducting 
qnantum interference device (SQUID) to the internal hyperfine states of a trapped ultracold atomic 
Fermi gas. Here, the SQUID, or flux qubit, acts as an effective magnetic impurity that induces spin- 
flip scattering near the Fermi energies of the trapped gas. Although the ultracold gas and SQUID 
are at vastly different temperatures, the formation of a strongly correlated Kondo state between the 
two systems is found when the gas is cooled below the Kondo temperature; this temperature regime 
is within current experimental limits. Furthermore, the momentum distribution of the trapped 
fermions is calculated. We find that it clearly contains an experimental signature of this correlated 
state and the associated Kondo screening length. In addition to probing Kondo physics, the con¬ 
trollability of this system can be used to systematically explore the relaxation and equilibration of 
a strongly correlated system that has been initially prepared in a selected nonequilibrium state. 


In the mid-1960s Jun Kondo gave the first detailed 
theoretical explanation of the anomalously large low- 
temperature resistivity that was observed in some metals 
[1]. This effect, which now bears his name, occurs when 
conduction electrons scatter off localized magnetic impu¬ 
rities with internal spin degrees of freedom. He showed 
that this exchange interaction leads to a breakdown of 
perturbation theory below an energy scale: the so-called 
Kondo temperature. This phenomenon was later under¬ 
stood in detail within Wilson’s renormalization group 
[2, 3]. In the renormalization sense, below the Kondo 
temperature an initial arbitrarily weak exchange coupling 
between the impurity and conduction electrons flows to 
the strong coupling regime. In the antiferromagnetic 
case, this flow ultimately produces a quasi-bound impu¬ 
rity state at the Fermi energy. It is the formation of this 
bound state that gives raise to the breakdown of pertur¬ 
bation theory and the enhancement of resistive scattering 

[4]. 

The advent of heterogeneous semiconductor devices 
brought about a dramatic resurgence of interest in the 
Kondo effect in the late-90s [5, 6]. In such systems, ar¬ 
tificial magnetic impurities atoms are created, manipu¬ 
lated, and coupled to Fermi gas leads. This enabled the 
controlled study of single and multiple Kondo impurities 
[7] in parameter regimes previously unexplored, includ¬ 
ing: magnetic field dependence, nonequilibrium effects, 
and non-Fermi liquid ground states [8]. More recently, 
scanning tunneling microscopes (STM) [9, 10] and spin- 
polarized STMs [11] have imaged the Kondo effect pro¬ 
duced by single magnetic adatoms, chains [12], and cor¬ 
rals [13]. While not as intrinsically tunable as quantum 
dots, coupling the impurity to an exotic bath state, such 
as a spin-imbalanced or superconductive state, is more 
easily realizable in adatoms systems [14]. 

Here, we propose a hybrid Kondo system whose ev¬ 
ery aspect is essentially fully controllable and tunable. 
The dimensionally, lattice structure, spin imbalance, in¬ 
teractions, and superconducting state of the Fermi gas 



FIG. 1. Schematic of the proposed composite Kondo system. 
The low-energy Hamiltonian of the SQUID can be modeled 
by an effective spin-Y 2 system with Zeeman splitting A. The 
two spin states of the Fermi gas correspond to different to¬ 
tal atomic angular momentum states, with a hyperfine energy 
splitting i?hf. The magnetic field B{r) produced by the clock¬ 
wise jo) and counter-clockwise jO) supercurrent states of the 
SQUID induces an exchange interaction JeS between the flux 
qubit and Fermi gas. When on resonance, A « Em, this in¬ 
teraction efhciently produces spin flips between the impurity 
and bath states via the Rabi cycle. 


can all be externally manipulated. The spin impurity 
itself is also tunable. Its effective Zeeman splitting can 
be dynamically controlled, as well as its position on the 
Bloch sphere. Furthermore, the gas-impurity exchange 
interaction can be tuned from weak to strong coupling. 
The proposed system is formed by coupling a supercon¬ 
ducting quantum interference device (SQUID), or flux 
qubit, to a two-component ultracold atomic Fermi gas. 
The low-energy Hilbert space of the SQUID is described 
by an effective Zeeman split spin-Y 2 Hamiltonian. The 
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two effective spin states are actually symmetric and an¬ 
tisymmetric combinations of left and right circulating 
macroscopic superconducting currents. These currents 
produce a magnetic field and, when on resonance, can 
induce magnetic dipole transitions between two selected 
hyperfine spin states of the trapped gas, see Fig. 1. Be¬ 
low we show that the combined system can be described 
by a Kondo-like Hamiltonian. Similar setups involving 
an atomic Bose-Einstein condensate have been recently 
put forth and realized experimentally [15-17]. 

A SQUID is a superconducting ring interrupted by 
a thin tunneling barrier. An applied magnetic flux in¬ 
duces clockwise | 0 ) or counterclockwise | 0 ) screening 
supercurrents in response. The current states are not 
energy eigenstates, but instead are eigenstates of the 
current operator / = la^ with eigenvalues ±/, corre¬ 
sponding to the total current and direction of each state. 
The two lowest energy eigenstates of the SQUID are 
IDs = 2-VDlO) + 10)) and |t)s = 2-VD|0) - |0)). 
In this basis the Hamiltonian of the flux qubit can be 
written as {h, = 1) iJs = typically A ~ 1-IOGHz 
[18]. Here the energy level spacing A will remain fixed, 
but in practice can be dynamically controlled on sub¬ 
nanosecond times scales [19, 20]. The macroscopic mag¬ 
netic field B{r) generated by the circulating currents 
can be determined from the Biot-Savart law, given the 
squid’s geometry and current /. In the current basis 
the magnetic field operator is B{r) = B{r)az. 

The fermion-fermion interaction of the two-component 
trapped gas can be externally tuned (through a Fesh- 
bach resonance) from strongly to weakly interacting. For 
simplicity we’ll assume they are in the weakly interact¬ 
ing regime and above any superfluid transition temper¬ 
ature [21]. The two experimentally prepared hyperfine 
states are atomic states of total internal angular mo¬ 
mentum F = J + I with z-projection mp and cor¬ 
responding internal energies Wo- that we label with a 
pseudo-spin- 1/2 index a =tD- For example, in the 
two states could be |t) = [A = 7j2,mp = —712) and 
ID = \F = 9/2,mp = —7/2), which have a hyperfine 
splitting of Ehfs = « 1.3 GHz. 

An interaction between the SQUID and the atomic 
gas is induced by the magnetic field produced by the 
flux qubit. This coupling can induce magnetic-dipole 
driven spin flips between the two hyperfine spin states 
and the two spin states of the SQUID. Neglecting cou¬ 
pling to other internal states, this interaction is described 
by U = -E, 7 ,a'/dT-^U’’)Ma.<T'^a'(r)( 8 ).B(r), where /x 
is the Lande g-factor weighted sum of each internal an¬ 
gular momentum contribution to the total magnetic mo¬ 
ment of the atom. In the weak-field limit /jl « gpfj,^F, 
where /ig is the Bohr magneton. Decomposing the elec¬ 
tron field operators in terms of the mode operators of 
the trapped gas 'I'CT^(r) = J2n^n\i')cnl and changing 
to the energy eigenstate basis of the isolated flux qubit 
(It) and It)), the total Hamiltonian can then be expressed 


as 
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where = en + cua- are the single-particle fermion en¬ 
ergies and with = 

—gpf^^J dr (j)’^{r)B^{r)(j)^,{r) is the effective exchange 
coupling. Equation (1) is the Hamiltonian of a Kondo- 
like system with a highly anisotropic exchange interac¬ 
tion and an applied Zeeman field. Typically, a Zeeman 
field tends to suppress the Kondo effect. Here, the effec¬ 
tive Zeeman splitting A of the impurity states is essential 
to its formation in this system, as the spin flip scattering 
is primarily driven by the Rabi process. During a Rabi 
cycle, the effects of the interaction are maximized when 
the two systems are on resonance, i.e., when Uhfs = A, 
and highly suppressed when far-off resonance. Although 
the Hamiltonian, Eq. (I), describes a Kondo-like system, 
it is far from thermal equilibrium, as the SQUID’s typi¬ 
cal operating temperature is on the order of a millikelvin, 
while the Fermi gas can be in the nanokelvin regime. As 
we will show in the following, the Kondo regime, or tem¬ 
perature, is set by the temperature of the gas. 

The nonequilibrium contour-ordered bath Green’s 
function is „,^,(z, z') = -i{TcCr,Az)cl>,^^{z')), 

where C = C+ U G_ is the Keldysh time contour with 
upper G+ and lower G_ branches [22]. The expecta¬ 
tion value is with respect to the non-interacting density 
matrix, while the time-evolution of the operators is with 
respect to the full Hamiltonian [23] . The contour Green’s 
function can also be expressed in matrix form as 


,{t,t') 


Gn(T,n'a-'{t,t') Gna,n'\ 
Gna,n'<7'i^^^') Gna,n'cr'(t,t') j 


where G[t, t) is the standard time-ordered Green’s func¬ 
tion, G{t,t) the anti-time-ordered one, and G^{t,t') are 
the lesser and greater correlation functions. Furthermore, 
to enable standard diagrammatic techniques we fermion- 
ize the impurity spin operators [24, 25] [26]. 

Although the Kondo effect ultimately leads to a break¬ 
down of perturbation theory, perturbative results remain 
quantitatively correct down to energy scales on the or¬ 
der of the Kondo Temperature Tk. To this end, we sum 
a subset of the leading-order logarithmically divergent 
diagrams, maintaining particle-hole symmetry, for the 
fermion self-energy [27]. The approximate self-energy is 
found in the so-called T-matrix approximation, or lad¬ 
der series, in both the particle-particle and particle-hole 
channels, see Fig. 2. The rung of each ladder series de¬ 
scribes the scattering of a bath and impurity fermion 
or hole. Within the steady-state approximation and 
Fourier transformed, the Keldysh matrix components of 
the rung of the T-matrix in the particle-particle channel 
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FIG. 2. The diagrams included in the particle-particle and 
particle-hole symmetric ladder approximation, or T-matrix, 
to the self-energy of the Fermi gas. The solid line represent 
the fermion propagator of the gas, while the dashed line is the 
projected pseudo-fermion propagator of the impurity. The 
spin and state dependent interaction between the two is de¬ 
noted by the wavy line. The subtraction of the last diagram 
in the self-energy is to avoid over counting the leading order 
term, which can be expressed either channel. 


where and in the particle- 

hole channel ^ (^ + 

)Gi^ (uj'), where is the non-interacting Green’s 

function of the bath and is the projected non¬ 
interacting pseudo-fermion Green’s function. 

The fermion self-energy in this approximation is ex¬ 
plicitly given by 
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where the T-matrices satisfy the following integral equa¬ 
tions in Keldsyh space 
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Within this approximation the the Kondo temperature 
will be taken as the temperature at which the T-matrices 
diverge; this signals the breakdown of perturbation the¬ 
ory. To make further progress we neglect the state de¬ 
pendence of the exchange interaction; —>■ / « 

where J is taken to be the value 
of the coupling for states near the Fermi surface and V is 
the volume of the trapped gas. Additionally, for now we 
set the squid’s temperature to be much larger than it’s 
level spacing. This essentially makes the two impurity 
states degenerate. One can then show that the above 
perturbation series diverges when the temperature of the 
gas feuTg approaches /cbIk ep exp [—1/(2| J|po)]) where 
po is the single-particle density of states per spin at the 
Fermi energy ep- Estimating the exchange coupling J 
for a two-dimensional gas with with Fermi wavelength 
Ap « 1 pm and Fermi energy ep « 1 pK that is placed 
100 pm above a 1 pm radius SQUID carrying 1 mA of cur¬ 
rent gives /cbTk O.Olep, which is within current exper¬ 
imental limits. Although the inclusion of the momentum 
dependence can change the functional relationship be¬ 


tween the coupling and Kondo energy scale, this estimate 
should still be qualitatively correct for weak coupling. 
One remaining uncertainty is the effect of decoherence on 
the formation of the Kondo state [28]. A detailed treat¬ 
ment of this will be left for future investigation, a simple 
estimate can be made by noting that the energy relax¬ 
ation time of a current flux qubit Ti ^ 1 ps and the time 
required for the Kondo singlet to form Ik ~ (feB^k)”^ 
[29] are of the same order. Thus, it should be possible 
for correlations to form before the SQUID decoheres. 

The Kondo effect can manifest itself as a large in¬ 
crease in the scattering rate, or inverse lifetime, of the 
bath fermions near the Fermi energy as the tempera¬ 
ture approaches the Kondo regime. The scattering rate 
is proportional to the imaginary part of the retarded 
self-energy which can be obtained from the lessor 
and greater self-energies. Figure 3 shows the imaginary 
part of the retarded self-energy for the spin-j, atomic 
fermions in various temperature regimes of the gas and 
the SQUID [30]. As the temperature of the gas is low¬ 
ered into the Kondo regime, the lifetime of the quasipar¬ 
ticles near the Fermi surface dramatically decreases. In 
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FIG. 3. The imaginary part of the retarded self-energy, which 
is proportional to the scattering rate, for the spin-J, atomic 
fermions is shown for various temperature regimes with (3s 
and /3g being the inverse temperatures of the SQUID and 
atomic gas respectively. The top panel shows that as the tem¬ 
perature of the gas crosses the Kondo temperature a Kondo 
resonance forms near the Fermi energy ep = p. The Kondo 
temperature for the system is £f/Ik ~ 150. The bottom panel 
shows how the scattering rate depends on the temperature 
of the SQUID when the Fermi gas is well below the Kondo 
temperature. When A/3s ^ 1 the spin states of the SQUID 
are essentially degenerate. As the temperature of the SQUID 
is lowered this gives raise to an effective magnetization of the 
impurity. 


a typical condensed matter system this enhancement of 
the scattering rate would express itself as an increased 
resistance. Measuring such transport properties in ul¬ 
tracold atomic systems is very challenging, but the mo¬ 
mentum distribution of the gas can be easily obtained 
by free expansion [21]. We find the momentum dis¬ 
tribution also contains a signature of the Kondo effect 
and the so-called Kondo screening length [31]. Roughly, 
this length scale sets the spatial extent of the singlet 
state around the impurity. To date, this length has 
never been experimentally observed. The momentum 
distribution of the atomic gas is given in terms of the 
spectral function; na{k) = J k,uj)nF(Lu — 

where Acr,a'{k,k,uj) = —2ImG^®^^,^,(w) and np{uj) is 
the Fermi distribution of the gas. To leading order in the 
self-energy this can be written as 

---^ /da;Im[S];,*;*^(cu)]nF(a; - ^ 

( 2 ) 

Figure 4 shows the momentum distribution above and be- 
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FIG. 4. The change in the momentum distribution of the 
atomic gas induced by Kondo correlations for a system far 
above the Kondo temperature (dashed line) and well below 
(solid line). In real space these correlations have a spatial 
extent given by the so-called Kondo screening length = 
vfI3k, where uf is the Fermi velocity and /3 k = {kBTK)~^. 
The associated momentum scale [fepl = 27r/^K defines the 
thin shell of momentum states near the Fermi surface which 
take part in the screening, i.e., those states with k — kpdikK. 


low the Kondo temperature. A qualitative understanding 
of these results can be had as follows: Deep in the Kondo 
regime the imaginary part of the self-energy is strongly 
peaked near the Fermi energy. The width of this peak 
is of order the inverse Kondo temperature Thus, 

the integral in Eq. (2) will only vary significantly when 
k « fcp ± fcx, where |feK| = 27r/^K and = vfPk is 
the Kondo screening length. Therefore, only momentum 
states within this shell of the Fermi surface participate 
in forming the Kondo cloud. 

We have shown that a fully controllable hybrid Kondo 
system can be formed by coupling a flux qubit to an ultra¬ 
cold atomic Fermi gas. The high degree of tunability of 
both the bath and the impurity offers the prospect of ex¬ 
ploring the Kondo effect in regimes never before realized. 
Critical to this—we find that the Kondo temperature for 
this system is within current experimental limits. The 
detection of the strongly correlated state between the 
impurity and bath could be achieved by measuring the 
momentum distribution of the gas or possibly by other 
methods such as rf-spectroscopy of the SQUID or atomic 
gas. This composite system also allows for a myriad of 
other interesting possibilities. For example, this system 
could be used to probe the thermalization of quantum 
systems [32]. The flux qubit can be initially prepared 
in an arbitrary state on the Bloch sphere then brought 
into resonance with the Fermi gas. This would allow the 
systematic study of the flow of the initial nonequilibrium 
state toward the correlated Kondo thermal state. 
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